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TECHNICAL MEMORANDUM X-64503 


GRAVITY GRADIENT MEASUREMENTS WITH A LASER 
ABSOLUTE GRAVIMETER 

SUMMARY 


This study presents the derivation of a technique for measuring gravity gradient, as well as 
absolute gravity, by modification of an existing laser absolute gravimeter. The technique has the addi- 
tional advantage of identifying a point in space corresponding to the obtained value of g. 

The method of error analysis used to predict the effect on performance resulting from known sources 
of measurement error is outlined. This analysis, as employed by the authors, indicates that expected 
accuracy is equal to that obtainable with the present state of the art. 

Brief consideration is given to the modifications which would be necessary to make the existing 
gravimeter compatible with the proposed technique. 


INTRODUCTION 


This study proposes an experiment designed to measure the rate of change of the acceleration 
caused by gravity (hereafter called “gradient” or “K”) through the employment of a laser interferometry 
technique. Such a technique is presently being utilized in a device for the measurement of the absolute 
acceleration caused by gravity. In essence, extensions of the theory employed in such a gravimeter are 
put forth, and modifications to its configuration which would render it amenable to gradient measure- 
ment are suggested. 

BACKGROUND MATERIAL 


Subsequent discussion relies heavily upon an article concerning laser absolute gravimetry 
which appeared in Laser Focus [1] . 

A laser system for measuring absolute gravity is shown in Figure 1 . The instrument is 
in the form of a Michelson interferometer in which mirror Ml is fixed and mirror 
M2 is a freely-falling body in the gravitational field g. Interference fringes are 
detected by the photodetector. Let x be the distance fallen along the gravitational 
field vector by the mirror M2 . . . ; then, at the photodetector, interference maxima 
are observed for: 


2x = 0,X,2X,3X, . . NX 




Figure L Basic relationship of primary components in laser interferometry 
system for measuring absolute gravity. 

The number of maxima, N, observed in a distance x is: 

2x = NX 


N = — x 

X 

Lambda is the wavelength of the illumination in a hard vacuum. If the light falling on 
M2 happens to make a small angle 6 with the gravity vector (remember that M2 
always falls parallel with the gravity vector), then the component of displacement 
parallel to g is x cos 0 , and the expression becomes: 

N = — x cos 0 
X 

Differentiating this expression twice with respect to time, we obtain a relation for g: 


-g = x = — - — N 
2 cos0 


[Note here that displacements downward along the gravity vector are considered 
positive and that the quantity g is inherently negative.] Subject to the initial con- 
ditions: N = 0, N = N 0 at t = 0 and integrating, one finds: 
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In the author’s experiment, time is determined by counting the number of oscillations 
of a standard oscillator of frequency f (hereafter called “the clock”). If n is the 
clock pulse number, then time is determined from: 


n = j fdt 
o 



a) 


Two electronic counters are used, one to record waves from the standard frequency 
oscillator and the other to record interference maxima from the photodetector. In 
this article, only integral values of n are considered. It should be noted that the 
higher the clock frequency f, the less time is associated with the integer n (see 
Equation 1). In an actual experiment, n is not really an integer, so the less time an 
integer represents, the less error it can introduce. The equation may now be written: 

n "(^)” s+ (t)“ < 2 > 

In Equation 2, the unknowns which must be experimentally computed are g and N 0 ’ 
To compute g, the quantity of interest, two equations are required to permit a simul- 
taneous solution. Hence, two measurements must be made for each experiment. Let 
us designate these measurements as A and B. Then, eliminating N 0 we find: 


g = 


_£A\/ N A n B' N B n A \ 
cos dj \ n A n B (n A + n B ) ) 


( 3 ) 


DISCUSSION OF INITIAL PROBLEM OF INTEREST 


This study was initiated by the problem of determining the point along the path of M2 at 
which the computed g value is applicable. This question arises from the fact that g is in reality not 
constant over the interval of the drop path. Its value is continuously increasing as M2 approaches the 
end of the drop. This can be seen from the Universal Law of Gravitation: 


F = -G Mm' 


m'g = -G 


mm 


g - - G 


m 

i 2 
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where m is the mass of the earth; m' is the mass of M2; G is the gravitational constant; r is the 
distance from earth’s center of mass; and g is the acceleration due to the earth. Thus, as M2 falls, 
it decreases the value of r which in turn causes an increase in the value of g. Differentiating the last 
expression above with respect to distance yields: 


_dg_ = + 2G -S- s K (1) 

dr r 2 

This quantity i s the gravity gradient, K. It constitutes the basis for the free-air correction 
dr 

utilized in geodesy to reduce gravity measurements to the surface of the accepted reference ellipsoid 
or spheroid. The universally accepted normal value of K utilized for this purpose is 3.08 X Iff 6 gal/cm. 1 

It can be seen that g is assumed to be constant throughout the drop path of M2 [1] . From 
the preceding discussion, it is evident that the value, obtained is merely representative of those values 
existing along the path. Then, to what point along the path does this representative value correspond? 

The effective length of the drop path utilized by the device discussed in Reference 1 is taken to 
be 16.593492 cm. Therefore, the difference between the g values at the beginning and the end of the 
path is of the order of 5 X 10“* 5 gal. Thus, if g measurements are made to seven figures or less, the 
value obtained is essentially representative of the entire path since the gradient effects do not enter until 
the eighth figure. 

Because advancing the state of the art will necessitate extending the accuracy of g measure- 
ments beyond the seventh figure, the problem of gradient effects is definitely of interest. So the 
problem of determining to what point along the path the measured g value corresponds becomes 
real. 


DERIVING A NEW EQUATION OF MOTION 


Since the whole problem centers around the assumption that g is constant over intervals of 
vertical displacement, the derivation of a new equation of motion in which g varies with distance is in 
order. In this derivation the gravity gradient, K, will be the constant of motion. First, however, it 
seems advisable to determine whether K can be considered constant over sufficiently large intervals 
that the problem incurred by considering g to be constant does not occur again. If we refer to 
equation (1) and differentiate once with respect to distance, we find 

JK = . go JS- = - 1.445 X 10“ 14 gal/cm 2 
dr r 


1. 1 gal = 1 cm/sec 2 . 
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where 


Gm = 3.986032 X 10 20 cm 3 /sec 2 


and 


r = 6.378165 X 10 s cm 


Thus, over a 1-m interval, K is seen to change only on the order of about 0.5 ppm. 

Taking K to be constant and displacement downward along the gravity vector to be negative, 
we begin our derivation. 

By definition 



Then 



dt 

dx 


= x 


1 

X 


JX 

X 
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We have previously defined as K; thus 

dx 

K = 4 (2) 

x 

The solution proposed for equation (2) is 
X = A + B e’^ + C 


Differentiating this expression to obtain its first three derivatives yields: 


x = 

x = A K + B K e"^ 
x = A K 3 / 2 - B K 3 / 2 e'^ 1 


Substituting the first and third derivatives into equation (2) shows that the solution satisfies the 
differential equation. Subjecting the solution and its derivatives to the initial conditions x 0 = 0, 

Kx 0 = x 0 , and x 0 = go at t = 0, the resulting system can be solved for the constants A, B, and C. 
This results in the following equation of motion: 


x 



sinh \/K(t) 




( 3 ) 


A similar result was derived by A. H. Cook [2] . 


EXPERIMENTAL SOLUTION FOR QUANTITIES OF INTEREST 


In equation (3), the unknowns to be computed from experimental data are x 0 ,g 0 , and K. 
Three equations are needed to afford simultaneous solution. Thus, three measurements must be made 
for each experiment instead of two as employed in the experiment discussed in Reference 1 . Designating 
these measurements as A, B, and C, one then arrives at the following system of equations: 
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VKx 


OA 


A K 
, _ ^ x OA 



gQA 

K 


*C 


K 

VKx 


cosh VK(t A ) - 1 


+ ^QA 
K 


£ sinh -\/K (t a + tg) - sinh \/K(t A ) 

— |^sinh \/K(t A + tg + t c ) - sinh \/K(t A + tg) 

+ cosh \/K(t A + tg + t c ) - cosh x/K(t A + tg) 

K L 


cosh -\/K(t A + tg) - cosh \/K(f ^) 


( 4 ) 


(The complete mathematical derivation of equations (4) may be found in Appendix A.) As the nature 
of the system makes any explicit solution for the unknowns impossible, one is forced to employ 
iterative techniques for solution. Solving the first equation for xq A and substituting this value into 

the second and third equations results in: 


‘ g 0A coshVK(t A ) ^ g 0A 


A K K 

j^sinh v^(t A + tg) - sinh %/K(t A ) 

sinh \/K(t A ) 


+ 


sinh \/K(t A + tg + t(-t) - sinh \/K(t A + tg) 
+ £ cosh \/K(t A + tg + t^O - cosh \/K(t A + tg)J 


g 0A cosh VK(t A ) 


*e 


K 


gQA 

K 


sinh \/K(t A ) 


gQA 

K 


cosh \/K(t A + tg) - cosi 


ihVK(t A ) 


Solving both expressions for gg A and subtracting one from the other, one obtains 


0 = K 


sinh \/K(t A ) 


A' 


*c- 


sinh Vf(t A ) 


C' 


T — V ’ - J*' 

B'- A' tanh y/K t D - C 'tanh \/K t ^ 
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where 


A' = sinh \/K(t^ + tg) - sinh >/K(t A ) 

B' = cosh VK(t^ + tg) - cosh \/K(t A ) 

C' = sinh VK(t A + tg + 1^) - sinh \/K(t A + tg) 

D = cosh \/K(t A + tg + t^) - cosh \/K (t a + tg) 

Thus, there results a function 0(K, t A , tg, Xq, x a , Xg, x^) = 0. This immediately lends itself to an 

iterative solution for K when values for the six measured parameters are substituted. The value of 
gQA is readily calculated from the results of that computation. 

Thus, one has obtained the values for two fundamental quantities, the absolute acceleration 
caused by gravity and the gravity gradient. In addition to this multiplicity of function, the value of 
g()A * s assi 8 ne d *o a particular point in space, xq A , the starting point of the first measured interval 

of distance. 


PERFORMANCE ANALYSIS 


To predict the performance of this experiment, and thus determine its feasibility, one must 
estimate the effect of measurement errors. This analysis must be performed on equations (4). The 
technique employed, the detailed explanation of which is found in Appendix B, was suggested by 
Charles Dalton of Aero-Astrodynamics Laboratory, Marshall Space Flight Center. Equations (4) are 
first put in the following form: 


F = 


Vk> 


OA. 


K 


sinhVK(t A ) + coshVK(t A )-l -x A = 0 

_ _ K 


G = 


Vk> 


0A 


K 


sinh %/K(t A + tg) - sinh \/K(t A ) 


g 0A 

K 


jcosh \/K(t A + tg) 
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H = 


y/Kx 

K 


OA 


|^sinh \/Kj(ty^ + tg + - sinh \/K(t^ + tg) 


jcosh + t B + *(9 " cos ^ V^^A + ' *C = 0 


where F,G, and H are all functions of K, Xq^, gg^, x^, Xg, x^., t^, tg, and t^. One then obtains 
the partial derivatives of the three dependent variables K, Xq^, and gg^ with respect to any particular 
independent variable by setting up a system of linear equations: 

Letting ue |x A , x B , t A , tg, t c | > 


3F = 

= _9F 

3K . 

3F 

9x 0A 

+ 9F 

9 80 A 

3u 

3K 

3n 

9x 0A 

3u 

9 80 A 

3u 

3G = 

s 9G 

3K . 

3G 

9x oa 

+ 3G 

9 §oa 

3u 

3K 

3u 

9x 0A 

3u 

9 §oa 

3u 

3H 3 

= _9H 

3K . 

3H 

9x oa 

, 3H 

9 §oa 

3u 

3K 

3u 

9x 0A 

3u 

a S0A 

3u 


(Expressions for the coefficients and constant terms involved may be found in Appendix C.) Thus, in 
matrix notation one has the following equation: 


3F 

3F 

3F 

3K 

9x 0A 

9 S0A 

3G 

3G 

3G 

3K 

^OA 

9 §oa 

3H 

3H 

3H 

3K 

9x 0A 

9 §oa 


3K 


3F 

3u 


3u 

9x oa 


3G 

3u 


3u 

9 §oa 


3H 

3u 


3u 
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Solving one such system for each of the independent variables will yield values for the partials of each 
dependent variable with respect to each independent variable. To obtain the probable error (P.E.) in 
any one of the three dependent variables, one needs only to substitute the values of the partials into the 
following equation: 


Letting ve j K, x 0A , g QA j 


fey m - fey w -fcyKT * &) w < 5 > 


The above technique may be used to predict performance, provided one has the values for all 
nine variables involved. As it is not possible to obtain the time values associated with each drop segment 
without constructing an apparatus to do so, solution of equations (4) for the dependent variables 
K, Xq^, and gg^ was out of the question for generating values to be used in the error analysis. 

Therefore, the values for x^, Xg, x^, K, Xq^, and g g^ were substituted into equations (4) and solved 
for times t^, tg, and t^. The displacements x^, Xg, and x^ were each taken to be -10 cm (roughly 
the length of segments used in the apparatus discussed in Reference 1). The value of gg^ was taken 

as -979.6395 gal, the value on a base plate in a nearby laboratory. The normal value 3.086 X 10~ 6 gal/ cm 
was used for K; and for purposes of simplification, xq^ was taken to be zero. Substituting these 

values into equations (4) and solving for the time values, one then has values for all nine variables suitable 
for purposes of error analysis. 


For the probable errors associated with the measured parameters x^, Xg, x^., t^, tg, and t^, 

we suppose the errors imposed by the limitations of the measuring devices employed in the apparatus 
discussed in Reference 1. For the distance measurements, one is, therefore, dependent upon the output 
of a gas laser of wavelength 6329.9147 A, which has a long-term stability of lCf 8 , and upon a fringe 
counter which is accurate to one-hundredth of a fringe. The time measurements will be derived from 
electronic counting of wave numbers (accurate to within one count) from a 100-Mc oscillator having a 
stability of 1(T 8 . 


With values for all nine variables and the necessary information concerning probable errors in 
the measured parameters, one has the data for computing the partials and P.E.’s contained on the right 
side of equation (5). The accuracy of these computations is admittedly limited by that of the extended 
precision mode of the IBM 1 130 computer that was employed. The data obtained from the analysis 
indicate that the performance will be well within required limits. 


RESULTS OF PERFORMANCE ANALYSIS 


It appears that a six-figure repeatability can be expected in the gg^ values obtained, which is 

the state of the art as attained by H. Preston-Thomas [3] and D. R. Tate [4] . In addition, five-figure 
repeatability in gradient values can be expected. (Sources known to the writers indicate that this would 
constitute a one-figure advance over present state of the art.) 
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To check the effect of varying values of K, gQ^, and x (x = = x^) in generating 

time values and the resulting effect on the performance analysis, trials were made for every possible 
combination of the following terms: 

K = 2.8 X 10" 6 gal/cm , 2.9 X 10 -6 gal/cm , 3.0 X 10~ 6 gal/cm , 3.1 X 10" 6 gal/cm 
gQA = -977 gal , -978 gal , -979 gal , -980 gal 
x = -5 cm, -10 cm, -15 cm, -20 cm 

The results for the probable error inherent to gQ^ values remained the same as in the first analysis. 

However, it was found that in some worst case conditions, the probable error inherent to K reduced 
the number of repeatable figures in gradient by one, to four repeatable figures. 

CONCLUSIONS 


This study indicates that the apparatus discussed in Reference 1 can be modified to make both 
gravity and gravity gradient measurements to state-of-the-art accuracies. The modification entailed 
would consist of little more than the incorporation of an additional segment of drop path and an 
additional counter. It would, however, also be of concern to determine the point at which the effective 
drop path begins, as the g value measured would correspond to that point in space. Location of this 
point within a millimeter seems to be possible by mechanical means. 
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APPENDIX A 


DERIVATION OF EQUATIONS (4) 

By definition 



x 


The proposed solution of the above differential equation is 
x = A e A ^ t + B + C 
Differentiating, 

x = A K ln eV Kt . B K m e'^ 1 
x = AKe^ + BKe^ 1 
x = A K 3/2 e^ 1 - B K 3/2 e"^ 1 


(A-l) 


(A-2) 
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Initial conditions at t = 0 are x 0 = 0, Kx 0 = x 0 , and x 0 - go- Imposing them on equations (A-l) and 
(A-2), one has 

0 = A + B + C 

Kx 0 = A K 3/2 - B K 3/2 

go = AK+BK 

Setting up the augmented matrix for tills system and reducing, one finds: 


1 

1 

0 

1 0 


1 

i 

1 

1 

1 0 
I 


1 1 

1 

1 

1 0 

k 3/2 

-k 3/2 

0 

! Kx 0 


0 

-2K 3/2 

o 

] 

J K^0-K 1/2 g 0 

I 


0 

1 

0 

1 rr 1/2 ' 

‘ K x 0 - g 0 








[ 2K 

K 

K 

0 

1 g0 


K 

K 

0 

1 

1 go 


1 

1 

0 

j go 

1 ~K 

] 


l 

i 


1 


1 k 1/2 x 0 + go 

1 1 1 l 0 


1 1 1 1 0 


1 o 0 I - 

1 

i 


1 

1 


1 2K 

1 

t ,1/2 * 

i K X 0 - g 0 


1 

1 IS 1/2 ' 

1 X X 0 - g 0 


i 

1 r/1/2 * 

K x 0 - go 

cs 

o 

o 

i 

cs 

o 

o 


° 1 ° ' - 2K 

1 

1 

J go 


1 

1 go 


1 

1 go 

0 0 - 1 | — 


L° 0 1 : j 


1 

o 

o 

1 


which yields: 

VKx 0 + g 0 

A = 

2K 
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B = - 


\/Kxi 


0 - go 


2K 


C 


go 

K 


The original system of equations thus becomes 



= / \ZKx 0 + go j y/Kt _ /ySxo_jo_\ .^/Kt 


^ Kxq + \/Kgo ^ e% /K t+ ^Kx 0 -VKgo j e _ > yK t 


They may be simplified to the following equations: 


x 


V^Xo / e VKt. e -VKt\ + / e VKt + e -VKt\ 80 
2K \ J 2K \ / K 



sinh(\/Kt) + — cosh(\/Kt) - 
K 


go 

K 


i - L'rtt t e-VK^ t — f • eVKt 

2K \ / 2K \ 

. _ \/Kgo 

= x 0 cosh(\/Kt) + sinh(\/Kt) 

K. 
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X = ( e V*. e -V*t) + _| ^ e VKt + e -VKt^ 

= y/Kx 0 sinh(\/Kt) + g 0 cosh(\/Kt) 

X = + e^Kt) + ^ ^ e v^t . e VKt j 

= Kx 0 cosh(\/Kt) + \/Kgo sinh(\/Kt) 

The equation of motion of an object during the first of three contiguous intervals would thus 
be 

X A = V ^° A sinh(V^t A ) + cosh(VKt A ) - 

K K K 

The equation of motion during the second of the three contiguous intervals is 
X B = N/ ^^ 0B sinh(^t B ) + cosh(VKt B ) - ^ 

where 

X 0B = X A = X 0A cosh (VKt A ) + sinh(VKt A ) 

K. 

SOB = X A = V^OA sinh(\/Kt A ) + g QA cosh(VKt A ) 
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Thus, the equation for Xg is 


Xg = sinh(VKtg) j^x 0A cosh(VKt A )+ ~^ A sinh(VKt A ) j 

+ -1- cosh(VKtg) ^ \/Kxq A sinh(s/Kt a ) + g QA cosh(\/Kt A ) J 
- -i- x/KxqA sinh (VKt A ) + g 0A cosh(\/Kt A ) J 

_ OA sinh(\/Ktg) cosh(\/Kt A ) + sinh(\/Ktg) sinh(\/Kt A ) 

K K 

+ cosh(\/Ktg) sinh(\/Kt A ) + cosh(\/Ktg) cosh(\/Kt A ) 

K K 

. V^OA sinh(N /Kt A) - cosh(VKt A ) 

Iv K 


VKx, 


ecx r 

sinh(\/Ktg) cosh(VKt A ) + cosh(VKtg) sinh(s/Kt A ) - sinh(\/Kt A ) 


g oA r 

K L 


sinh(v^tg) sinh(\/Kt A ) + cosh(\/Ktg) cosh(\/Kt A ) - cosh(\/Kt A ) 


Tliis results in the following equation: 
n/KxqA 


X B = 


K 


sinh(\/Kt A + \/Ktg) - sinh(\/Kt A ) 


+ [ cosh(\/Kt A + \/Ktg) - cosh(%/Kt A ) J 


The equation of motion during the last of the three intervals is 


*C 


- ^SL sinh(VKt c ) + cosh(\/Kt c ) - M 
K K K. 
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where 


x 0 c = x B = x 0B cosh(VKt B ) + sinh(VKt B ) 

K. 


gQC ” X B ” x/KxqB s i n K\/Kt B ) + gQB C0S h(VKt B ) 


Xqb = X A = x 0A cosh(v/Kt A ) + V ^ 0A sinh(v/Kt A ) 

K. 


SOB = *A = sinh(s/Kt A ) + g 0A cosh(v/Kt A ) 


Xqc =[ x OA C0S KVKt A ) + sinh(VKt A )l cosh(v/Kt B ) 


+ VK . 

+ sin 

K 


inh(>/Kt B ) ^x/Kxqa sinh(VKt A ) + g 0A cosh(\/Kt A ) 


gQC = VK sinh(VKtg) Xq A cosh(\/Kt A ) + s inh(v/Kt A ) 

K 


+ cosh(\/Ktg) >/Kx qa sinh(>/Kt a ) + g QA cosh(VKt A ) J 


Reducing, one finds that 


x OC = X 0A cos K\/Kt A ) cosh(\/Kt B ) + s inh(\/Kt A ) cosh(\/Kt B ) 

K 


xq A sinh(\/Kt A ) sinh(\/Kt B ) + s inh(\/Kt B ) cosh(\/Kt A ) 

K 
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X 0C - X 0A 


cosh(\/Kt^) cosh(\/Kt j|) + sinh(\/Ktg) sinh(\/Kt A ) 


VKg, 


OA 


K 


sinh(\/Kt^) cosh(\/Ktg) + sinh(\/Ktg) cosh(\/Kt A ) 


|+ ^gQA 


x QA cosh(VKt A + VKtg) + sinh(VKt A + VKtg) 

K 


and 


gQc = VKxqA sinh (VKt B ) cosh(VKt A ) + g QA sinh(VKt B ) sinh(VKt A ) 

+ \/Kxq A cosh(\/Kt B ) sinh(\/Kt A ) + gQ A cosh(\/Kt B ) cosh(\/Kt A ) 
= VKxqA j^ sin h(\/Ktg) cosh(\/Kt A ) + cosh(VKtg) sinh(VKt A ) J 
+ gQ A j^sinh(\/Ktg) sinh(VKt A ) + cosh(\/Kt B ) cosh(\/Kt A )J 
- %/Kxqa sinh(VKt A + \/Kt B ) + g QA cosh(VETt A + \/Kt B ) 
Substituting these values in the original equation for x^, one finds that 


Xq = sinh(\/KtQ) 


,, ^ 0 A 


x 0A cosh(VKt A + VKt B )+ -'^ a - sinh(VKt A +VKt B ) 

K 


K 


CO; 


sh(VKt c ) 


\/Kxq A sinh(\/Kt A + VKtg) + g 0A cosh(\/Kt A + VKtg) 


1 

K 


VKx qa sinh(VKt A + VKtg) + g 0A cosh(\/Kt A + VKtg) 
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X£ = sinh(\/Kt^) cosh(\/Kt^ + \/Ktg) + ^2^ sinh(\/Kt^) sinh(\/Kt^ + \/Ktg) 

K ^ 


K 


K 


+ cosh(\/Kt^) sinh(\/Kt A + \/Ktg) + cosh(\/Kt^) cosh(\/Kt A + \/Ktg) 

K 


. s inh(VKt A + \/Ktg) - C osh(VKt A + VKtg) 

K K 


Further simplifying, one finds that 


*C = 


vfe, 


OA 


K 


g OA 

K 


sinh(\/Kt^) cosh(\/Kt^ + \/Ktg) + cosh(\/K t^) sinh(\/Kt^ + VKtg) 


- sinh(\/Kt^ + VKtg) 


sinh(\/Kt^) sinh(\/Kt^ + \/Ktg) + cosh(\/Kt^) cosh(\/Kt^ + VKtg) 


■ cosh(\/Kt^ + VKtg) 


VKx, 


OA 


K 


g OA 

K 


sinh(\/Kt^ + VKtg + V®^) - sinh(\/Kt^ + VKtg) 


cosh(VKt^ + VKtg + VKt(V - cosh(\/Kt^ + \/Ktg) 


Thus, the resulting system is 

Vk^qa 


Xa = 


Xd = 


K 

n/K>L 


sinh VK(t A ) 


gQA 

K 


cosh VK(t A ) - 1 


OA 


K 


inh \/K(t A + tg) - sinh \/K(t A ) J + jcosh \/K(t A + tg) 

-coshVK(t A )] 
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xq = ' £ sin h VK(t A + tg + tj-.) - sinh VK(t A + tg) 

+ — ^ cosh VK(t A + tg + 1^) - cosh %/K(t A + tg) J 
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APPENDIX B 


MEMORANDUM CONCERNING ERROR ANALYSIS TECHNIQUE 


The following is the text of a memorandum by Charles C. Dalton on the subject “Error 
Propagation in Gravimetry Experiments.” This information was provided in response to a specific 
request by the authors of this paper for technical assistance. 

The mathematical aspects of the problem, as it now is, is that three parameters to be 
calculated, say x, y, and z for initial velocity, gravity value at t = 0 and gravity 
gradient, are related to six measured parameters, say a, b, c, d, e, and f by three 
simultaneous equations of some considerable complication, say 


F(x, y, z, a, b, c, d, e, f) = 0 

(1) 

G(x, y, z, a, b, c, d, e, f) = 0 

( 2 ) 

H(x, y, z, a, b,c, d, e, f) = Q 

( 3 ) 


You want to solve the equations for x, y, and z and to compute the random error 
of each propagated from the measured parameters. 

An explicit solution for x, y, and z might be complicated and difficult to find or 
even not possible. But the given implicit equations can be programmed for a numeri- 
cal solution to whatever accuracy is required; e.g., geometrically each of the equations 
can be considered a surface in the three dimensions x,y, and z. The solution, which 
for physical reasons is known to exist, must be one of the possibly several points of 
intersection of the three surfaces. Calculating the propagated error does not require 
the solution to be explicit. 

The six measured parameters are three measurements of time and three measurements 
of distance in terms of wavelengths of the laser light. So far as I can tell. . . the random 
errors of the six measured parameters can be considered effectively both statistically 
independent and all normally distributed. . . . [Some further thought should be given 
to this point.] 

The formulation for the propagation of error parallels that of the total differential; i.e., 

dx = — da + — db + — dc + — dd + . . . 
da 3b 3c 3d 
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where an acceptable approximation is to replace each of the differentials by corre- 
sponding small finite “delta” increments. But this presupposes that both the 
(arithmetic) magnitude and the algebraic sign of the increments of the independent 
parameters are known, whereas the random errors in specific measurements are 
unknown in both aspects. Therefore, with statistical independence, the different 
terms of the total differential are added like orthogonal vector components by 
squaring each component; and the finite “delta” increments can each be expressed as 
the same small fraction A of the corresponding standard deviation o, say Aa a , 
Aa^, etc. Then, by multiplying both sides of the equation by A*" 2 


o 2 = (3x/3a) 2 a a 2 + (3x/3b) 2 a^ 2 + (3x/3c) 2 o 2 + ... (4) 


and similarly for y and z. But without statistical independence there would also 
be other terms involving correlation factors. The result, rigorous only for linear 
functions, is an approximation which omits, usually, smaller terms with higher 
derivatives. 

In principle, the six partial derivatives in equation (4) could be evaluated by 
giving small increments alternately to the six independent measured parameters and 
noting the change to the solution for x, y, and z by the given simultaneous 
equations (1) through (3). Undoubtedly this would require special double precision 
programming and more computation than would be required by solving for each of 
the six measured parameters a, b, . . . a set of three linear equations relating the 
partial derivatives of the three dependent variables x, y, and z with respect to the 
particular dependent variable, e.g., a; i.e., by partial differentiation of equation (1), 


0X 

3a 


+ (9F/9y) + (9F/9z) 

3a 3a 3a 

_ 3F 
3x 


If the correctness of this manner of differentiation of an implicit function is not 
obvious, it may be helpful to rewrite the same result in symmetric form; i.e., 

(9F/9x) — + (9F/9y) + (3F/9z) ^ ^ = 0 . 

9a 9a 9a 9a 
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By differentiating equations (2) and (3) similarly and switching to subscripts to 
denote partial differentiation, say 

F„ for and x Q for 

x 9x a 3a 


one gets three linear equations for the partial derivatives of the three variables x, y, 
and z with respect to the measured parameter a; i.e., in matrix notation 



( 5 ) 


and similarly for each of the other five measured parameters b, c, etc., by replacing 

a. Note that the coefficient matrix is the same for all of the measured parameters a, 

b, etc. 



APPENDIX C 


EXPRESSIONS FOR TERMS INVOLVED IN ERROR ANALYSIS 


The terms involved in error analysis are as follows. 


9 F = 
9 K 


= 2 ^ g 0 A l A ' X 0 A sinhVK(t A ) + ^ t A - M C oshVK(t A ) + 


ff = gOAV^A sinhVK(t A + t B )-sinhVK(t A ) 

+ ~2K tA "K^ C ° Sh + tfi) ' C ° Sh %/K ' (t A ) 


2S_ \i 

2 K 


cosh >/K(t A + tg) + 


sinh \/Ki (t^ + tg) 


U = gQA^A + ■ X 0 A sinh VK(t A + t B + t c )- sinh VK(t A + t B ) 


(t A + ’ ~“T cosh A + *B + ' cosh A + 


+ -^- x 0A cosh VK(t A + t B + t c ) + > ^ 0A sinh\/K(t A + t B + t c ) 


9 F = ^JT 

3x 0A k 


9 G = 

9Xn a 


^ |^sinhVK(t A + t B ) - sinh V^t A )] 


_ 9 H_ _ \/E_ sinh Vi 5 (t A + t B + t c )- sinh \^(t A + t B ) 
,X 0 A K L 
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9F _ 1 

9g 0A K 

j^cosh \/K(t A ) - lj 


9G _ 1 

3gQA K 

jcosh \/K(t A + tg) - cosh \/K(t^)j 


9H _ 1 

9 §oa k 

j^cosh \/K(t^ + tg + t^) - cosh \/K(t A + tg) 


S’ SJ 

II 

9F =0 

9 X g 

3F =0 

9xc 

9G = 0 
9x a 

9G = j 
9x B 

9G = 0 
9x c 

o 

li 

pc < 
«> 

3H = o 
9x B 

9H _ j 

9x c 

9F _ • 
^ ° A 

cosh VK(t A )j + k ° A j^sinh\/K(t A ) 


9t^ = * 0A 
at A 

cosh \/K(t ^ + tg) - cosh VK(t^) + j^sinh \/K(t^ + tg) 



- sinh VK(t^) 

9H _ • 
at 0A 

3t A 

cosh \/K(t^ + tg + t^,) - cosh \/K(t^ + tg) 



+ sinh VK(t^ + tg + t^.) - sinh \/K(t^ + tg) 

K 
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9F _ 


at 


= o 


B 


_3G 

dtn 


= X, 


OA 


cosh \/K(t^ + tg) 


VKg, 


OA 


K 


sinh \/K(t A + tg 


9H _ • 

at B 0A 


cosh %/K(t A + tg + t^) - cosh \/K(t A + tg) 


VKg, 


OA 


K 


sinh \/lC(t A + tg + t^) - sinh \/K(t A + tg) 


_3F 

dtp 


= o 


dG _ 


atr 


= 0 


dH 

dtp 


= x. 


OA 


cosh \/K(t^ + tg + t^) 


VKg, 


OA 


K 


sinh \/K(t^ ■ 


tg + t c ) 
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